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ABSTRACT 

The spatial and spectral properties of three solar flare coronal X-ray loops are studied before, 
during and after the peak X-ray emission. Using observations from the Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI) , we deduce the temporal changes in emitting X-ray length, 
corpulence, volume, position, number density and thermal pressure. We observe a decrease in 
the loop length, width and volume before the X-ray peak, and an increasing number density and 
thermal pressure. After the X-ray peak, volume increases and loop corpulence grows due to an 
increasing width. The volume variations are more pronounced than the position variations, often 
known as magnetic line contraction. We believe this is the first dedicated study of the temporal 
evolution of X-ray loop lengths and widths. Collectively, the observations also show for the first 
time three temporal phases given by peaks in temperature. X-ray emission and thermal pressure, 
with minimum volume coinciding with the X-ray peak. Although the volume of the flaring plasma 
decreases before the peak in X-ray emission, the relationship between temperature and volume 
does not support simple compressive heating in a collapsing magnetic trap model. Within a low 
/? plasma, shrinking loop widths perpendicular to the guiding field can be explained by squeezing 
the magnetic field threading the region. Plasma heating leads to chromospheric evaporation and 
growing number density, producing increasing thermal pressure and decreasing loop lengths as 
electrons interact at shorter distances and we believe after the X-ray peak, the increasing loop 
corpulence. 

Subject headings: Sun: Flares - Sun: X-rays, gamma rays - Sun: corona - Sun: chromosphere 



1. Introduction 

During a solar flare, it is believed that magnetic 
energy in the corona is converted to other forms 
of energy; causing the bulk motion and heating of 
surrounding plasma and accelerating particles up 
to non-thermal energies. If the energy is released 
in a low density coronal region, then it is possi- 
ble that the majority of accelerated particles will 
leave the corona without losing energy collision- 
ally until they reach the much higher densities 
of the underlying chromosphere. Here the par- 
ticles will interact with the surrounding plasma 
and produce hard X-ray (HXR) bremsstrahlung 
emission. However, in most flares we also ob- 
serve coronal X-ray emission at or near to the top 



of a l oop (see iHolman et al.l 120111 : iKontar et al 
2011aL as recent reviews), although with an in- 



teresting exception ( Fleishman et al.ll201l[ ). Since 
bremsstrahlung emission grows with density, the 
low coronal densities can often make the obser- 
vation of coronal loop emission challenging, espe- 
cially if it is masked by brighter chromospheric 
X-ray emission (Krucker & Lin .2008,) . However, 
the observation of HXR coronal sources with in- 
struments such as the Ramaty Hi gh Energy Solar 



Spectrometer Imager (RHESSI) ( Lin et al.l 120021 ) 



are important since these X-ray loops are much 
closer to or may even contain the actual sites of 
energy release and primary electron acceleration 
(IXu et all [20081: [Kontar et al.l lioTlbl : iBian et al ' 
20121 : IGuo et all |2012|) . RHESSI creates X 
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ray images us ing Rotating Modu l ation Collima- 



tors (RMCs) (jHurford et alJ 120021 : ISchmahl et al 



20071 ). As the spacecraft spins, the grids placed in 



front of the detectors produce a modulation curve 
of the X-ray intensity. From this we can construct 
X-ray visibilities, which are equiva lent to the vis- 



ibilit ies used in radio astronomy (jSchmahl et al 



20071 ). In the last few years, an imaging technique 



known as visibility forward fitting (VIS FWD- 
FIT) has been successfully used to s tudy the sizes 



of chromosp heric footpoint sources |Kontar et al 



[2008. 201 0: Battagha fc Kontail2011l ) and inferred 
that the cross-sectional sizes of HXR sources de- 
crease with energy. This technique has the advan- 
tage over other imaging algorithms such as Clean 
and Pixon etc. because it provides us with rel- 
atively accurate position and size measurements 
from the moments of the X-ray distribution, plus 
an error for each value. However, X-ray visibility 
fitting also has the disadvantage of limiting the 
number of sources and their shapes that can be 
studied using this technique. Therefore, careful 
comparisons with other imaging algorithms are 
required when using VIS FWDFIT. 

Recently, a number of studies have used visibil- 
ity forward fitting to study t he siz es of coronal X- 



ray lo op sources. IXu et al.l (|2008l ) and lGuo et al 



(|2012n investigated how the lengths of X-ray coro- 
nal loops, that is, the direction parallel to the 
guiding magnetic field, varied with photon energy 
during the X-ray impulsive stage. For a num- 
ber of different fiares, they found that the loop 
lengths increased with energy with a form con- 
sistent with that of an initial extended accelera- 
tion region within the loop itself, plus an addi- 
tional length proportional to the photon energy 
squared. This is due to higher energy electrons 
moving further throu g h the plasma before inter- 



acting. iKontar et al.l (|2011bl ) studied the length 



and width changes of one X-ray coronal loop with 
photon energy. They found that both the length 
and width increased with photon energy, with the 
loop width increasing proportionally to the pho- 
ton energy. Increases in loop width, or corpulence, 
are more difficult to explain since the electrons are 
bound to the guiding fi eld and cross fi e ld tran sport 
should be neglig ible. Kontar et al. ( 2011b[ ) and 



Bian et al.l (|2011i ) inferred that the width increase 



could be due to the magnetic diffusion of field lines 
perpendicular to the direction of the field, caused 



by the presence of magnetic turbulence within the 
loop. Although, these observations indicate the 
usefulness of coronal loop spatial properties, the 
spatial evolution of loop sources has not yet been 
fully explored. Loop positions have been stud- 
ied before, d uring and after the imp ulsive phase 
of the flare. Forbes fc ActonI ( 1996[ l observed in 
soft X-rays the changing locations of post flare 
loops, which was interpreted as the decrease in 
height that open fleld lines undergo after they have 
reconnected to form closed loops. They studied 
two long duration events near the limb and found 
that field loop shrinkage did occur and matched 
the shrinkage predicted by a simple model of the 
reconnecting field but found ove rall the entire 
fiare l o op system grew with time. Sui &: Holman 
(l2003l).ISui et al.l (l2004l ). IVeronig et all (|2006l ) and 



Joshi et all ( 20091 ) all noted a decrease in the al- 



titude of loop top sources during the impulsive 
phase of the flare, until the peak X-ray emission 
and a n increase in a l titude after the impulsive 
phase. ISui fc HolmanI (|2003l) also found evidence 



for an above the loop top source and interpreted 
the situation as the formation of a reconnection 
current sheet between th e loop top sou r ce and the 



higher coronal source. IVeronig et al. (l2006l) 



terpre ted the decrease as a c ollapsing magnetic 
trap ( Somov fc Kosugi 1997 : Karlicky fc Kosugl 
20041 ). The contraction and expansion of the 



loop source has also been observed in other wave - 
lengths of EUV (iLiu et alJ '2009': J oshi et al.ll2009l) 
and radio (Li fc Can 2005 : iReznikova et al.l201o[) . 
Reznikova et al. ( 2010( ) found both changes in the 
radio l oop spa n and height with time. More re- 



cently, ICosain ( 2012 ) looked for evidence of col- 
lapsing flelds using Solar Dynamics Observatory 
(SDO/AIA and HMI) observations. The loops 
rose slowly and then moved into a collapse phase 
during the impulsive phase of the flare, where the 
loop tops contracted. Lower loops contracted ear- 
lier than higher loops and the loop contraction was 
interpreted as a reduction of magnetic energy as 
the system relaxe d to a state o f lower energy, i.e. 
relaxation theory (|Tavloiill974l ). 

In this paper, we study three flares with coro- 
nal X-ray loops to find how the emission lengths, 
widths and positions change with time at three 
energy ranges between 10-25 keV. From each flare 
spectra, we will also investigate how parameters 
such as emission measure and plasma temperature 
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Table 1 

Table showing the main parameters of each flare. 





GOES Class 


Date 


Obs. time 


Peak time (10 keV) 


Footpoints (30-40 keV) 


Flare 1 


M3.0 


23-August-2005 


14:22:00-14:40:00 


14:30:00 


14:36:00 onwards 


Flare 2 


M4.1 


14/15- April-2002 


23:58:00-00:20:00 


00:12:00 


00:05:00 onwards 


Flare 3 


M2.6 


21-May-2004 


23:42:00-23:58:00 


23:50:00 


23:42:00 onwards 



vary with time. Using a combination of imaging 
and spectroscopy parameters, we can infer how 
the X-ray loop volume, number density, thermal 
pressure and energy density vary during the time 
evolution of the flare, allowing us to build a fuller 
picture and propose some explanations describing 
the processes occurring within the coronal loops. 

2. Chosen events with coronal X-ray emis- 
sion 

The three events studied are: 23rd August 2005 
from 14:22:00 (Flare 1), 14th/15th April 2002 
from 23:58:00 (Flare 2) and 21st May 2004 from 
23:40:00 (Flare 3). All three flares share similar 
characteristics: GOES M-class flares with similar 
lightcurves, strong coronal X-ray loop top emis- 
sion and only relatively weak footpoint emission. 
Since the aim of our study is to examine how the 
properties of coronal loop emission vary with time, 
these events were chosen as they show loop emis- 
sion throughout the rise, peak and decay stages of 
X-ray emission and their s patial properties have 
been previou s ly stud ied by IXu et al.1 (|2008t ) and 
Kontar et al.l (|2011bl ). The coronal X-ray emis- 
sion during each flare appears up to ~ 25 keV and 
each source is a simple loop-like shape connect- 
ing HXR footpoints at 30-40 keV during certain 
time intervals. The main parameters of each flare: 
GOES class, date, observation time, peak time 
at 10 keV and the time of footpoint appearance 
are given in Table [TJ Length variations of each of 
these coronal loops wi th photon energy were stud- 
ied bv lXu et"aL ( 20081 ). Length and width changes 
with p hoton energy for Fla re 2 were also stud- 
ied by iKontar et al. ( 2011b( ). It should be noted 
that Flare 1 and Flare 3 show similar results as 
Flare 2 in iKontar et al.l (|2011bl ) , where both the 
length and corpulence increase with energy as ^ 
and ^ e respectively but this paper will concen- 



trate on the size, position and spectral parameter 
changes with time. The lightcurves for each event 
are shown in Figure [5] (top row) . Flares 1 and 3 
have similar lightcurves; a simple shape with one 
peak. The lightcurve for Flare 1 is shown in Fig- 
ure [S] (top row, left plot). Our study of this event 
begins at 14:22:00. At this time, emission from the 
10 to 20 keV energy bands are slowly rising and 
reach a peak at ^14:30:00. After this point, there 
is a gradual decrease in X-ray emission, which can 
be observed up until ~14:50:00, where RHESSI en- 
ters into a night phase. In the 20-40 keV band we 
see a series of peaks between our observation range 
of 14:22:00 and 14:40:00. The lightcurve for Flare 
3 is shown in Figure [5] (top, right plot). From the 
start of our study at 23:40:00, the X-ray emission 
is rising from 6-50 keV and peaks around 23:50:00. 
After this peak, the X-ray emission decreases. The 
lightcurve for Flare 2 is more complex and is shown 
in Figure [5] (top row, middle plot). During the ob- 
servational time there are two main peaks in the 
lightcurve at -00:03:00 and -00:12:00, possibly 
more peaks, which are easiest to see in the 20-40 
keV energy range. 

3. Spatial and spectral changes with time 

Each event was studied using the imaging al 
gorithms of Clean (Hogbom Il974 : iHurford et al 



2002), Pixon (|Pina fc Puetteilll993HMetcalf et al 



19961) and VI S FWDFIT (iHurford et all l2002t 



Schmahl et al.l 120071 ). Firstly, we studied each of 



our three events using only Clean and Pixon to 
confirm the loop shape of each coronal source and 
find the energy ranges over which a coronal source 
was present in each flare. Once we were confl- 
dent that our chosen events only had a simple loop 
shape, we studied each event using VIS FWD- 
FIT. This is important since the coronal source 
must have a loop-like shape so that VIS FWD- 
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Fig. 1.— 23-August-2005 (Flare 1). Clean image (green) at 14:28:00-14:30:00 and over-plotted VIS FWDFIT 
contours (50% maximum intensity) at four selected times for 10-12 keV (left), 12-15 keV (middle) and 15-20 
keV (right). The asterisks denote the loop position for each selected time. The cyan VIS FWDFIT contour 
matches the time of the Clean image. 




Fig. 2.— As Figure [U but for Flare 2, 14th April 2002. 




Fig. 3. — As Figure [T] but for Flare 3, 21st May 2004. Flare 3 uses the energy ranges of 14-16 keV (left), 
16-20 keV (middle) and 20-25 keV (right). 
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FIT can effectively fit a curved elliptical gaussian 
(loop) to the X-ray visibilities of the event and 
give realistic estimates with errors for the source 
parameters. VIS FWDFIT provides us with loop 
length FWHM (full width half maximum), loop 
width FWHM and the {x, y) centroid position of 
the loop. It should also be noted that VIS FWD- 
FIT gives the mean position of the loop shape, not 
the loop top centroid position, that is, the central 
position at the top of the loop. The coordinates for 
loop top position were obtained by extracting the 
coordinates of the central circular gaussian that 
VIS FWDFIT uses to create the final loop shape 
along with a set of other circular gaussians placed 
along the length of the loop. This is important 
as a loop that is very curved and approaching the 
shape of a ring will pull the shape centroid to- 
wards the ends of the loop, often masking small 
changes in loop position with time or energy. This 
was especially significant for the large loop shape 
of the 23rd August 2005 event (Flare 1). In or- 
der to study changes in time at a specific energy 
range, we imaged the coronal sources of Flares 1, 2 
and 3 over two minute intervals where possible and 
four minute intervals where the count rates were 
lower. The exact time bins used for each flare are 
shown in Figures [5] and |6l The energy ranges of 
10-12 keV, 12-15 keV and 15-20 keV were chosen 
for Flare 1 and Flare 2 and the energy ranges of 
14-16 keV, 16-20 keV and 20-25 keV were chosen 
for Flare 3. These energy ranges were chosen since 
VIS FWDFIT with a loop produced clear results 
with small errors at these energy ranges. Images 
for Flares 1, 2 and 3 are shown in Figures[Tl[2]and[3] 
respectively. Figure [1] shows that Flare 1 is a limb 
event while Figures [2] and [3] show that both Flare 
2 and Flare 3 are disk events (the main parame- 
ters of each flare are given in Table [J) . Figures 
[H [2] and [3] plot a background Clean image of the 
coronal source for each flare at one chosen time 
interval and over plot VIS FWDFIT contours for 
selected time intervals, one corresponding to the 
same time interval as the Clean image. Comparing 
the shape and size of the VIS FWDFIT contour 
with the Clean background image at the selected 
time interval for each flare shows good agreement 
between both algorithms. The Clean images here 
use a sigma beam width of 1.8. It should be noted 
that we also found the standard deviation from 
image cuts through the images of Clean, Pixon 



and VIS FWDFIT and it was found that a Clean 
sigma beam width of ~3.0 showed the best agree- 
ment with the standard deviations of Pixon and 
VIS FWDFIT. In terms of imaging, VIS FWD- 
FIT also agrees well with the Clean and Pixon 
algorithms at all other times not shown in this pa- 
per. From Figure [SI it should be noted that the 
results for Flare 3 are probably less reliable than 
those of Flare 1 and Flare 2. From the position 
of the footpoints in the Clean image, it appears 
as through the southern 'loop leg' is tucked un- 
derneath the observer's line of sight. This means 
that it is harder for VIS FWDFIT with a loop 
to fit it with a correctly shaped loop and usually 
fits it with a loop that is slightly too large or de- 
parts from a loop-shape. For each of the imaging 
time intervals for Flares 1, 2 and 3, energy spectra 
were created. The spectra of each flare at each 
time interval were fitted with a thermal compo- 
nent and a non-thermal component correspond- 
ing to thick target bremsstrahlung. The thermal 
fits provide us with the emission measure, EM , 
and the plasma temperature, T. We note all three 
flares show the emergence of footpoints in the 30- 
40 keV energy range. Spectra for Flares 1, 2 and 
3 are shown in Figure Slat three selected time bins 
corresponding to a rise, peak and decay stage of 
X-ray emission for each flare. 

Figure [5] shows the imaging parameters for 
Flare 1 (column 1), 2 (column 2) and 3 (column 3): 
loop width FWHM (row 2), loop length FWHM 
(row 3) and loop-top radial position (row 4) for 
each imaging energy band. Above these plots, the 
lightcurve is plotted for each of the imaging en- 
ergy bands to allow comparison with changes in 
the spatial parameters. The dashed lines drawn 
in each plot represent the time interval over which 
the peak X-ray emission occurs for our three en- 
ergy bands. For each flare, in general, peaks in 
the lightcurve represent changes of width, length 
and source position parameters with time. The 
last two rows of Figure [5] plot how the spectro- 
scopic parameters vary with time for each flare. 
In general, for all three events, the emission mea- 
sure rises throughout the observation time, either 
slowing or decreasing slightly during the last few 
minutes for each event. For each event, the plasma 
temperature peaks before the peak in X-ray emis- 
sion and the n slowly decreases aft e r this point, 
as noted by lAntiochos &: Sturroc^ (|l978l ). The 



5 



plasma temperature decreases much slower than if 
it were decreasing by conduction only, suggesting 
prolonged energy release at the later decay stages 
of the flare. 

3.1. Loop corpulence 

For all three flares, we see a general pattern 
for emitting loop width changes with time. Be- 
fore a peak in X-ray emission, the source width 
for each flare tends to decrease and after a peak 
in X-ray emission, the source width for each flare 
increases. It is also notable that the rates of both 
width expansion and contraction (width change in 
time) are approximately the same as can be seen 
in Figure [71 

Flare 1- For all energies plotted, the source 
width decreases until the peak X-ray emission at 
-14:30:00-14:32:00 UT and then increases after 
this point. The largest change occurs for the high- 
est energy of 15-20 keV. This falls from 14.5" to 
11.5" at 14:28:00-14:30:00 and then increases untfl 
14:40:00 where the width peaks at ~ 20", produc- 
ing a final larger loop corpulence than seen at the 
beginning of our observational time. 

Flare 2- For the energies of 10-12 keV and 12-15 
keV, the loop width decreases before the first peak 
in X-ray emission at —00:03:00. After this point 
the width increases before dropping at 00:12:00- 
00:13:00 where there is another X-ray peak in the 
lightcurve. After this peak, the width continues to 
grow again. Before 00:03:00, the width of the 10- 
12 keV source falls from around —10" to — 5" and 
after 00:13:00 rises from -6" to -13". The 15-20 
keV source also shows this pattern except there 
is a larger peak at 00:05:00 and then a more pro- 
nounced decrease in width until 00:12:00-00:13:00 
UT. 

Flare 3- Again, the change in loop width with 
time follows a similar pattern as Flares 1 and 2. 
For all three energies ranges considered, the loop 
width decreases from 23:42:00 to the peak in X- 
ray emission at 23:50:00. The loop width then 
increases after this time. The 14-16 keV and 16-20 
keV sources - faU from 7-8" at 23:42:00 to -5" at 
23:50:00 and then increase up to 8-9" at 23:58:00. 
The 20-25 keV source falls from 11" to 5". Note 
the missing data at the fourth and sixth time bins 
for the 20-25 keV energy range, as VIS FWDFIT 
was unable to fit successfully at these times for 



this energy range. 

3.2. Loop length 

As with loop width, each fiare shows general 
pattern for emitting loop length changes with 
time. Before a peak in X-ray emission, the source 
length for each flare tends to decrease and after a 
peak in X-ray emission, the source length either 
increases slightly or remains approximately con- 
stant within the errors. 

Flare 1 - For all energies plotted, there is 
a rapid decrease in source length until the peak 
X-ray emission at -14:30:00-14:32:00. After the 
peak, the source length remains approximately 
constant (within the error) . The smallest decrease 
in loop length before the X-ray peak emission oc- 
curs for the 10-12 keV source which falls from 54" 
to 38". The decrease in loop length grows with en- 
ergy and the highest drop in loop length occurs for 
the highest energies plotted at 15-20 keV, which 
fall from 74" to 38" before the peak in X-ray emis- 
sion. 

Flare 2 - The pattern for length changes with 
time are similar to that of the width changes. 
For the 10-12 keV source, the length before the 
00:03:00 peak in the lightcurve, falls from —26" 
to -20", rises to -24" at 00:11:00, falls to -20" 
at 00:15:00 and then increases to 21" at 00:17:00. 
Due to the multiple peaks, the loop length tends 
to increase after a peak in the lightcurve. 

Flare 3 - For all energies, the length drops 
rapidly between 23:42:00 and 23:50:00 from the 
range of 25" to 30" to —10" for all energies plotted. 
Again as with Flare 1, after the X-ray peak the 
length of the loop remains approximately steady 
until the final plotted time of 23:58:00 for all en- 
ergies. 

3.3. Loop radial position 

For all three flares, peaks in each X-ray 
lightcurve tend to denote times where the trends 
in loop radial position change. 

Flare 1- Since this is a limb event, the radial 
position can tell us whether the source is moving 
away or towards the limb. Before the X-ray emis- 
sion peak at 14:30:00-14:32:00, the source moves 
towards the limb, falling a distance of — 2" at 
10-12 keV, 12-15 keV and 15-20 keV. Then after 
the peak, the source moves away from the limb. 
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Fig. 4. — Spectra for Flares 1 (top row), 2 (middle row) and 3 (bottom row), at three chosen imaging time 
bins (during the X-ray rise, peak and decay stages). Black — data - background, Olive — background, red 
= thermal fit, green = thick target bremsstrahlung fit. The residuals are plotted below each plot. 
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By plotting the actual source positions, we also 
see that overall the entire loop structure moves 
northwards during the time interval of 14:22:00- 
14:40:00, moving in a U-shape as also shown for a 
number of flares in Shen et al. I (l2008l) . From our 
results for Flare 1 , we see that the changes in radial 
position with time are '^comparable to the width 
changes and smaller than the length changes. The 
largest change in position is only 2 — 3" while the 
width decreases by 3" and increases by 6". The 
length shows the largest change with a decrease of 
~ 17" or more before the peak in X-ray emission. 

Flare 2 - Overall, the radial positions of the 
10-12 keV, 12-15 keV and 15-20 keV sources in- 
crease with time. At points of peak X-ray emis- 
sion, that is, at -00:03:00 and -00:12:00, we see 
changes in the gradient. The slope steepens be- 
tween -00:03:00 and 00:12:00. Since this is a disk 
event, it is difficult to say whether there is a change 
in source altitude at these peaks (as for Flare 1). 
By plotting the source positions on the solar disk 
with time, we can see that before —00:03:00 the 
source moves in a north-easterly direction, it then 
moves southwestward until —00:12:00, before mov- 
ing northeastwards after this peak. At 10-12 keV, 
between 23:58:00 and 00:20:00 the source radial 
direction changes by 9.5". Therefore the overall 
change in position is larger than the individual 
changes in loop width and loop length for Flare 2. 

Flare 3 - For all energies, the radial distance 
falls with increasing time. There does not seem to 
be any significant difference in the radial distance 
trend after the X-ray emission peak at 23:50:00, 
apart from the steadier decrease in radial distance 
after this time. As with Flare 2, Flare 3 is a 
disk source and hence it is difficult to determine 
if there is any altitude change. Plotting the (a;, y) 
disk position shows us that the source ccntroid 
moves in a north-easterly direction until the peak 
at —23:50:00 and then begins to move in a south- 
easterly direction, once there is a decrease in the 
X-ray emission. For all three energies, the radial 
distance decreases by —8" between 23:40:00 and 
23:58:00. 

4. Changes in emitting volume and other 
inferred parameters 

From the loop width FWHM W and loop length 
FWHM L, we can infer the general changes in the 



emitting loop volume, V , for each flare over time, 
at each energy band. We assume that the vol- 
ume of the loop is given hy V = ttW'^L/A that is, 
assuming a cylindrical loop. The changes in emit- 
ting loop volume with time for Flares 1 , 2 and 3 are 
plotted in Figure|6] (second row). In general, for all 
three events, the source volume decreases before a 
peak in X-ray emission and increases after a peak 
in X-ray emission. The changes in plasma num- 
ber density, thermal pressure and thermal energy 
density can all be calculated using combinations of 
loop volume, emission measure and plasma tem- 
perature. The plasma number density, n, can be 
obtained via n = y'EM/V , the pressure, P, from 
P = nksT, where fcs is the Boltzmann constant 
and finally the energy density, E = 3nkBT. The 
variation of these quantities with time for Flares 
1, 2 and 3 are shown in Figure |6] (third, fourth and 
fifth rows respectively). 

Flare 1 - As expected from the width and length 
results, the loop volume falls between 14:22:00 and 
the peak in the X-ray lightcurve at —14:30:00 and 
then rises after this time for all three energies. For 
all three energy bands, the decrease and increase 
in loop volume occurs at roughly the same rate. 
The largest decrease is for the highest 15-20 keV 
band, faUing from — 4.7 x 10^"^ cm^ to — 1.4 x 10^'^ 
cm'^ at 14:29:00 and then rising after this time to 

— 4.7 X 10^^ cm^^ at the last observational time. 
The 10-12 keV band falls from - 2.4 x 10^'^ cm^ 
to - 1.3 X lO^^^ cm^ at 14:31:00 and then rises 
back to — 2.5 x 10^*" cm'^ at the final observa- 
tional time. The number density, thermal pres- 
sure and energy density for all three energy bands 
tend to follow the same pattern, rising to a peak 
at some time after the peak X-ray emission and 
then slowly decreasing. For the 10-12 keV band, 
the number density rises from 1.5 x 10"'^'^ cm~^ at 
14:22:00 to 6.5 x 10^° cm'^ at 14:35:00. It then 
falls to - 5.6 X 10^" cm-3 at 14:38:00. The 12-15 
keV and 15-20 keV bands follow similar patterns, 
peaking at —14:33:00. The pressure rises from 

- 40 g/[cm s^] at 14:22:00 and reaches 170 g/[cm 
s^] at 14:31:00, the time where the X-ray emission 
peaks. After this time the pressure remains ap- 
proximately constant at — 160 — 170 g/[cm s^] un- 
til the last observation time where it falls to — 140 
g/[cm s^]. The 12-15 keV and 15-20 keV energy 
bands peak at —14:33:00. The energy density is 
just the pressure multiplied by 3 and hence it fol- 
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Fig. 5.— (left) 23-August-2005, (middle) 14-April-2002 and (right) 21-May-2004. 1st row - lightcurves (see 
graph legend for energies), 2nd row - loop width vs. time, 3rd row - loop length vs. time, 4th row - radial 
position vs. time, 5th row - emission measure vs. time and 6th row - plasma temperature vs. time. Dashed 
lines = peak X-ray emission. 
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lows the same pattern as pressure throughout the 
observed duration of the flare. The energy density 
peaks at 500 ergs cm~'^ in the 12-15 keV band at 
~14:33:00, after the peak X-ray emission. 

Flare 2 - For the 10-12 keV and 12-15 keV 
sources, the loop volume falls until the first peak 
at 00:03:00 and then increases until it reaches 
00:12:00, drops at 00:13:00 and then increases 
again after this time. As for Flare 1, the num- 
ber density, thermal pressure and energy den- 
sity all follow the same pattern for each energy 
band, rising to a time at or just after the peak 
in X-ray emission and then decreasing after this 
point. The highest number density, thermal pres- 
sure and energy density for each energy band oc- 
curs at 00:12:00 to 00:14:00, just after the peak in 
X-ray at 00:11:00-00:13:00. The number density 
peaks at around 18 x 10^" cm~'^ in the 12-15 keV 
band, while the thermal pressure and energy den- 
sity peak at ~ 500 g/[cm s^] and 1500 ergs cm~^ 
respectively at this time and energy. 

Flare 3 - The loop volume again falls before the 
peak in X-ray emission at 23:50:00 and then rises 
again after this time, for all three energies. For 
this flare again, the number density, thermal pres- 
sure and energy density all peak just after the peak 
in X-ray emission. At 23:50:00-23:52:00 the num- 
ber density peaks at 23 x 10^° cm"'^, the thermal 
pressure peaks at 750 g/[cm s^] and the energy 
density at ~ 2300 ergs cm~'^, in the 16-20 keV 
band. 

5. Summary and discussion 

Using visibility forward fitting, we have for the 
first time produced a dedicated study of how the 
spatial and spectral properties of three coronal 
X-ray loops changed during the temporal evolu- 
tion of each flare. All X-ray loops exhibited sim- 
ilar changes in both their spatial and spectral 
properties and hence the results indicate that a 
common process is occurring for all three events; 
the emitting flaring loop volume is decreasing be- 
fore the peak in X-rays and increases after the 
peak in X-rays is reached. Before the peak X- 
ray emission, the emitting lengths and widths of 
each coronal loop decreased with time, indicating 
that the X-ray emitting region of the loop volume 
was contracting (a reduction in loop corpulence 
and length) as the X-ray emission from the region 



grew. After the X-ray peak, the loop corpulence 
increased at approximately the same rate as dur- 
ing contraction stage. For Flares 1 and 3, with 
one peak in their lightcurves, once the minimum 
X-ray loop length was reached during the X-ray 
peak it remained approximately constant, at least 
within the errors of our results. 

Sirnilar to prev ious studies /Anti ochos fc Sturrock 
19781: ISiding fc Spiccr 1980: Gunkler et al.lll984t 



McTiernan et al.l 1199^ spectroscopy for each 



event showed that the plasma temperature ini- 
tially grew but began to decrease before the peaks 
in X-ray emission and emission measure. The 
emission measure for each flare generally grew 
with dips observed during the flnal observational 
times of Flares 2 and 3. At the same time, the 
number density, thermal pressure and energy den- 
sity of the plasma also increased as the X-ray emis- 
sion grew. The plasma temperature decreased 
much slower than by thermal conduction only, 
even during the X-ray decay phase, suggesting at 
later stages additional energ y release is required 
(e.g. iKolomahski et al.l 120111 ) to explain the long 
lasting X-ray loop emission. For the limb event 
(Flare 1), we do see a decrease in loop altitude be- 
fore the peak in X-ray emission of roughly 2" ~ 1.4 
Mm, which is comparable in magnitude to the 
decrease in loop width but overall, the largest 
changes occurred for the X-ray loop volume, not 
the X-ray loop position. Decreases in loop altitude 
before the peak in X-ray emission have been well 
noted before and arc often referred to as coronal 
implosion, loop contraction or flcld line shrinking. 
Sui et al. ( 2004h also studied the loop position 



changes of Flare 2. From their observations, they 
did conclude that the loop was indeed decreasing 
in altitude before the first peak in the lightcurve at 
00:04:00. They concluded that the loop decreased 
by ^ 2" over a 4 minute period, which is consistent 
with our radial distance results for Flare 2. How- 
ever, no altitude decrease before the larger peak in 
X-ray emission at 00:12:00 was noted, where again 



we see decreases in loop length and loop width. 

It has been suggested that a decrease in loop al 
titude may be an indication of collapsing magnetic 
trap acceleration/heating (ITakakara & Kai 1966 




Somov & Kosugil 


1997 


; IPriest & ForbesI 12002 


Karlickv & KosuEil2004 


:IVeronig et al.ll2006:lGradv et al. 


2012). Although this mav be the case for other 



events, we do not believe that the collapsing mag- 
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Fig. 6. — (left column) 23rd August 2005, (middle column) 14th April 2002 and (right column) 21st May 
2004. 1st row - Lightcurves for each event, 2nd row - Volume changes with time, 3rd row - Number density 
changes with time, 4th row - Thermal pressure changes with time and 5th row - thermal energy density 
changes with time. 
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P.= 



027 + 0.002 kms"' 



016 + 0.002 krris"' 



14:24 14:28 1 4: J2 1 4:.36 14:40 
Start Time (23-Aug-05 14:22:00) 

Fig. 7.— For Flare 1 (23rd August 2005), 
lightcurve (top row), dW/dt (2nd row), dL/dt 
(third row) and dr/dt = v (bottom row). Both 
dW/dt and v can be fitted using a straight fine 
indicating constant accelerations aw and a„ and 
hence constant forces. 



netic trap is the prime solution to our observa- 
tions, firstly due to overah larger loop volume 
changes over relati vely small position changes. 
Veronig et al. 1 (l2006l) observed the coronal loop of 



a GOES X-class flare and found downward veloc- 
ities as large as ^ 14 km/s in 10-15 keV band 
and ~ 29 km/s in the 15-20 keV band prior to 
the peak in X-ray emission. This is not what we 
observe for these events. For Flare 1, we only ob- 
serve an average downward velocity in the 10-20 
keV band of ~ 4 km/s, which is generally com- 
parable to the decrease in loop corpulence during 
this period. Average changes in loop radial posi- 
tion for Flares 2 and 3 are also ~ 4 km/s. More 
convincingly, in a collapsing magnetic trap model, 
simple compressive heating during the contraction 



stage would imply that NT g<: 1/A (jMaetzler et al 



[l978; Enislic 198l|, where N is the number of par- 
ticles in the region, T is the plasma temperature 
and A is the cross-sectional area of the region, 
given by wW"^ , where W is the loop width. Figure 
|9] plots log NT against log A for each flare and 
shows this not to be the case. Straight lines fits to 
both the contraction and expansion phases show 
that the gradients are either greater or less than 
-1. 

Observations by iKontar et al.l (|2011bl ) showed 
how the loop corpulence of the 14/15th April 2002 
event (Flare 2) increased with energy during two 
times corresponding to the first rise and first de- 
cay stages of our observations and suggested that 
the presence of magnetic turbulence (the diffusion 
of field lines) was the cause of the energetic width 
increases. Hence, the suggestion that magnetic 
turbulence in the region can account for the ener- 
getic changes in loop corpulence, may also be able 
to ac count for the extra e nergy in the loop , since 
both lKontar et all (|2011bl) and lBian etaU (|201l[ ) 
found the energy density of magnetic fluctuations 
to be significant and could be comparable to that 
of the flaring plasma and higher than the energy 
density of non-thermal particles. Flares 1 and 3 
also show similar length and width increases with 
energy at a single time range, also suggesting the 
presence of magnetic turbulence. 

Ma ny observations ( Li fc GanI 2005t Liu et al 



12009; 'Joshi eLalJ 120091: iReznikova et al.l l2010t 



.Gosain 2012i) of loop height and length changes 
have been explained in terms of a reduction in 
magnetic pressure. Usually the reduction in mag- 
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Fig. 8.— SOHO EIT 1951 images for Flare 1 at the times of 14:21:12 and 14:34:51, corresponding to the 
times of rise and peak in X-ray emission. RHESSI 10-12 keV X-ray contours are over plotted in pink. Note 
the lack of bright EUV emission at the rise stage of this flare. 



netic pressure is referred to as Taylor relaxation 
but this only refers to a special c ase where th e 
resulting field is linear force free (j Taylor! Il974l) . 
The reduction in magnetic pressure could also 
account for the reduction of loop corpulence 
or cross-sectiona l area, as shown in simulations 
by IJanse fc Lowl (2003) and hence the observed 



trends of number density and pressure. iLiu et al 



(|2009[ ) studied coronal implosion of one coronal 
source. They explained the reduction in height of 
a coronal source before the peak in X-ray emission 
in terms of Taylor relaxation. They also suggested 
that this type of event can only occur if the coro- 
nal loop is already filled with hot, dense plasma 
before the onset of new event, that is from a pre- 
vious event in that region. The fact we do not 
see EUV emission during the rise phase of Flare 
1 seems to co rrelate wi t h the observations and 
suggestions of iLiu et all (|2009[ ). Figure |8] shows 
SOHO (Solar and Heliospheric Observatory) EIT 
(Extreme Ultraviolet Imaging Telescope) 195A 
images at the times of 14:21:12 and 14:34:51 for 
Flare 1. X-ray emission contours at 10-12 keV for 
14:22:00-14:26:00 and 14:34:00-14:36:00 are also 
over plotted. From Figure 13 we can see that dur- 



ing the rise phase there is no bright 195A EUV 
emission emanating from the loop, only 10-20 keV 
X-ray emission. After the peak in X-ray emission, 
we do see EUV emission from the loop. We see 
an overall increase in the number of particles, N, 
within the loop region throughout the duration 
of all three events. Due to low coronal densities, 
chromospheric evaporation probably accounts for 
this increasing iV, initially driven by thermal con- 
duction and possibly at later times by electrons 
reaching the chromosphere, when we see weak 
footpoint emission and EUV emission from the 
loop. 



5.1. Three temporal phases and suggested 
explanations for the observations 

For each flare. Figure [10] replots the plasma 
temperatures. X-ray emissions, emitting loop 
widths and thermal pressures but now at a single 
energy band of 10-20 keV (14-25 keV for Flare 3). 
From Figures [5] and [B] and now more clearly in 
Figure [TOl we can see that collectively, the obser- 
vations of each fiare display three distinct phases 
and each of these phases will form the basis of 
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Fig. 9. — Plots of logA^T against logl/A for Flares 1 (left), 2 (middle) and 3 (right). The star represents 
the start time of each event while the triangle represents the peaks in X-ray emission for each event. The red 
and blue dashed lines represent straight line fits during the compressive (red) and expansive (blue) phases. 
For Flare 2, only the first compressive and final expansion phase have been fitted. For simple compressive 
heating, we would expect the gradient of each line to be 7 = —1. 



our suggested explanation. Each phase is repre- 
sented by a shaded orange bar. During Phase 1, 
there is a peak in plasma temperature and during 
Phase 2, a peak in X-ray emission. At Phase 2 
we see the smallest emitting loop width, length 
and hence volume. Finally, during Phase 3, the 
thermal pressure of the region peaks. For Flare 
1 (left column), each phase is well separated and 
can be clearly seen. The pattern can also be seen 
for Flares 2 (middle column) and 3 (right column) 
but each phase is not as clearly defined as for 
Flare 1. During Flare 3, each phase occurs over 
much shorter time intervals and therefore each of 
the phases overlap slightly. During Flare 2, there 
are multiple peaks, which along with the shorter 
timescales for each process, makes each individual 
phase harder to see. However, the overall pat- 
tern is observed for all three fiares. We can see 
from Figure [10] that each phase can only be eas- 
ily seen for slower events. The quicker the event, 
the harder it is to distinguish between each of the 
three phases as each phase overlaps in time. Many 
solar eruptive events may exhibit a similar pattern 
but each phase may occur at timescales too fast 
to be currently observed. 

In order to understand our observations, we 
believe it is crucial to understand the decreasing 
X-ray widths before the X-ray peak. Length vari- 
ations along the field can always be explained in 
terms of changing number density, but in order to 
explain the corpulence variations, where electrons 



are tied to the magnetic field, the most plausible 
explanation is the movement or the diffusion of 
magnetic field lines. We have already noted the 
energetic loop width increases at a given time for 
each flare, suggesting the presence of turbulence 
in the region, but overall in time the width at a 
given energy band shrinks until the X-ray peak 
is reached, implying the presence of both turbu- 
lence and the shrinking of the cross-sectional area 
of the field at this stage of the flare. Therefore, 
we suggest the plasma within our emitting X- 
ray loop region is tied to the magnetic field lines 
and hence the contraction of the emitting loop 
cross-section or corpulence during the rise phase 
is ultimately due to the contraction of the cross- 
sectional area of field lines that thread the region 
or possibly the expansion of field lines above the 
region. It is also sensible to assume that the re- 
gion we call the loop in X-rays actually consists 
of multiple coronal loops that cannot be resolved 
in X-rays using RHESSI. Though we can only see 
what is happening to the region of the loop that 
is emitting X-rays, we will assume the entire loop 
region is contracting since the plasma is tied to 
the held lines. Hopefully, from future observa- 
tions we will find similar, newer coronal events 
with complimentary AIA/SDO data where it may 
be possible to see if the region consists of mul- 
tiple coronal loops. From our observations, the 
reason for this cross-sectional width shrinking is 
unknown and beyond the scope of this paper, but 
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Fig. 10. — Observations of plasma temperature, X-ray emission, loop width and thermal pressure are re- 
plotted together for Flares 1 (left), 2 (middle) and 3 (right) at one energy band of 10-20 keV (14-25 keV 
for Flare 3). The orange bars represent three phases: 1. a peak in plasma temperature, 2. a peak in X-ray 
emission, generally coinciding with the smallest loop width and 3. a peak in thermal pressure. 
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we can speculate for further study that it may be 
due to a reduction in magnetic pressure within 
the X-ray loop region, as often suggested to de- 
scribe height decreases (coronal i mplosion or loop 
contracti on) disc ussed earlier ( Li fc Gan _ 2005t 
Liu et al.„ 2009: Joshi et al.ll2009t iReznikova et al 



20101 : lGosainll20li r" The magnetic pressure de- 



creases as the magnetic field relaxes from a non- 
potential state. Although it is reasonable to as- 
sume that during a solar eruptive event such as a 
flare, the field will reside in a non-potential state. 
X-ray observations of these events do provide 
evidence for the non-potentiality, again through 
the inference of magnetic turbulence and hence 
a non-parallel field component within the loop 
region. The non-potential state and decreasing 
magnetic pressure could be due to reconnection 
abov e the region or maybe even along the loop it- 
self (IVlahos et al" 2004t Gordoyskvv fc Browning 



2OIII l2012t iGordovskvv et al.ll20l'^ . One such 
model is caused by a kink instability as shown 
in rec e nt sim u lations by Gordovskvv fc Brownine 



(|201lL 12012); IGordovskvv et all (|2012D . In this 



model, the energy of twisted loops is released by 
reconnection inside the loop and transferred to 
plasma heating and particle accelerations. 

During phase 1, the process causing the con- 
traction of loop width/cross-section is also proba- 
bly responsible for the comparable decrease in loop 
altitude, or at least the change in loop position ob- 
served for Flares 2 and 3 sitting on the solar disk. 
Plots otdW/dt (width contraction/expansion) and 
dr/dt — V (centroid velocity) for Flare 1 plotted 
in Figure [7] show that each parameter can be fit- 
ted using a straight line during both the contrac- 
tion and expansion phases. dL/dt (length contrac- 
tion/expansion) is also plotted in Figure [7] but its 
trend cannot be described by a straight line fit 
suggesting the changes in length are caused by a 
different process, itself a consequence of the pro- 
cess producing the decreasing loop corpulence. It 
is also sensible to assume that this process is re- 
sponsible for the increased plasma temperature 
and the acceleration of electrons. The increasing 
temperature of the region means that energy will 
be thermally conducted towards the lower levels of 
the solar atmosphere causing gentle chromospheric 
evaporation of the denser coronal and chromo- 
spheric layers below. The observations show the 
plasma temperature peaking relatively early be- 



fore the peak in X-ray emission and then slowly 
decreasing; slower than by thermal conduction, 
implying that energy is still being supplied to the 
loop plasma. This is most likely via the conver- 
sion of magnetic energy. Chromospheric evapora- 
tion drives plasma into the region producing the 
increasing number density and hence thermal pres- 
sure, along with the shrinking corpulence at this 
phase. The increasing number density is responsi- 
ble for the rapid non-linear decreasing X-ray loop 
length, since electrons accelerated within the re- 
gion will travel shorter distances before interact- 
ing. 

Phase 2 During Phase 2, after the peak in 
plasma temperature, the loop corpulence stops 
shrinking. Number density and thermal pres- 
sure within the loop still increases due to chromo- 
spheric evaporation and the length of the emitting 
region also reaches its lowest point. The loop cor- 
pulence may stop shrinking because the process 
causing the shrinking ceases or it may be due to 
the balancing of forces within the region. For ex- 
ample, if the reduction of loop corpulence is due 
to the reduction of magnetic pressure in the loop 
then at this phase, the growing thermal pressure 
may finally be high enough to balance the reduc- 
tion in B pressure. 

Phase 3 During the final phase, the thermal 
pressure continues to rise within the loop due to 
the increasing number density from chromospheric 
evaporation. We believe the growing thermal pres- 
sure in the region is now responsible for the ex- 
panding loop corpulence. This expansion, in turn, 
eventually halts the increasing number density and 
thermal pressure at a time after the peak in X- 
ray emission. After Phase 3, the loop corpulence 
continues to increase and we see slow decreases in 
both number density and thermal pressure. The 
X-ray emission continues to decrease during and 
after Phase 3 and the emitting loop length dur- 
ing this period remains approximately constant, 
equal to the minimum loop length in Phase 2, even 
with a decreasing number density. It is sensible 
to assume that the acceleration mechanism in the 
loop is slowing during this time. However, Flare 
2 is an exception to this trend as we see multiple 
events/X-ray peaks in the lightcurve. 

Finally, we would like to note that a stud y of 
the 23rd July 2002 flare bv lCaspi fc uA (|2010[ ) cal- 
culated temporal volume changes using the Clean 
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algorithm and assuming an elliptical geometry. 
Overall this flare shows a general trend consistent 
with our results; an overall decrease in volume be- 
fore the peak in X-ray emission and an overall in- 
crease in volume after the peak in X-ray emission. 
This flare also shows a peak in plasma tempera- 
ture before the peak in X-ray emission and a high 
number density (and hence thermal pressure) after 
the X-ray emission first peaks. 
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